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The elastic scattering differential cross section is calculated for proton scattering from ®He at 717 
MeV, using single scattering terms of the multiple scattering expansion of the total transition am- 
plitude (MST). We analyse the effects of different scattering frameworks, specifically the Factorized 
Impulse Approximation (FIA) and the Fixed Scatterer (adiabatic) Approximation (FSA) and the 
uncertainties associated with the use different structure models. 

PACS numbers: 24.10.-i, 24.10. Ht, 24.70. +s, 25.40. Cm 



INTRODUCTION 



The understanding of nuclei can only be fully achieved by studying the way they interact with other nuclei. This 
is particularly true for halo nuclei, due to their short lifetime and simplified energy spectrum. We consider here 
^ . the scattering of a structureless projectile by a target assumed to be composed of structureless subsystems (as is 
approximately the case for halo nuclei). 

It is the aim of a microscopic reaction theory to construct the total scattering amplitude in terms of well defined 
dynamical and structural quantities. For the scattering of a nucleon by a target with J\f particles, one has to solve a 
n = A/" + 1 many body scattering problem. For M —2 clusters of equal mass this has been done solving the Faddeev 
equations Pi and generaHzed for TV = 3 [1 Q ■ This many body scattering framework is however very complicated, 
' and does not handle clusters of different mass, and so alternative methods have been developed as for example the 
Continuum Discretized Coupled Channels (CDCC). In this approach a system of coupled equations needs to be solved 
with effective projectile-subsystem interactions. This method has been successfully appHed for more than two decades 
to the scattering of {N — 2)-body targets for a wide range of projectile masses and bombarding energies. The case 
of the scattering from a {N = 3)-body system is considerably more demanding, although some work is already in 
progress Q . Alternativelyin the high energy regime one can use a multiple scattering expansion of the total transition 
amplitude, MST [10,@,Bllli- 

When describing the scattering of stable from halo nuclei, it is crucial to model the halo many-body character 
of composite particles 0, |^, U- In particular, a three or four-body problem has to be solved when studying the 
scattering of a projectile by a target which is a bound state of two or three subsystems, as is the case for ^^Be and 
^He respectively. 

This problem can be conveniently addressed by the MST approach. In this formalism, the projectile-target transition 
amplitude is expanded in terms of off-shell transition amplitudes for projectile-subsystem scattering. Due to the 
complexity of the many-body operator, suitable approximations need to be made in order to express in a convenient 
way the overall scattering amplitude in terms of the scattering by each target subsystem. Under further suitable 
approximations each term can be written in terms of a product of a form factor and the transition amplitude for the 
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scattering for that subsystem. The MST method provides a clear and transparent interpretation of the scattering of a 
composite system in terms of the free scattering of its constituents and is numerically advantageous. This scattering 
framework is particularly useful at high energies and for loosely bound nuclei where the expansion is expected to 
converge quickly. 

In this paper, we use this framework to calculate the elastic scattering differential cross section for proton scattering 
on ^He at 717 MeV where new data has been obtained at higher transferred momentum than previously. These new 
momentum transfers, however, still essentially probe only that part of the few-body dynamics of the halo cluster which 
is constrained by the rms radius. Therefore the few-body treatment of the halo is appropriate in this energy and 
angular range, and we will see that calculated differential cross sections for various structure models largely reflects 
properties closely connected to rms radii. It is the aim of the present paper to clarify and test the approximations 
involved in the applying the MST multiple scattering expansion to find the predicted cross sections for proton-halo 
elastic scattering at these momentum transfers. Specifically, we examine (i) the impulse approximation, (ii) the 
factorization approximations and (iii) the single scattering approximation. 

II. MULTIPLE SCATTERING EXPANSION OF THE TOTAL TRANSITION AMPLITUDE 

We consider the scattering of a projectile (system 1) from a few-body target consisting of M sub-systems weakly 
bound to each other. We shall frequently refer to the composite system as the target although in practice an actual 
experiment may be carried out with the composite system as a projectile. The subsystems X, . . . are assumed to 
be stable and can be either composite nuclei or nucleons. The total transition amplitude for the scattering is 

T = V + VGT 

n 71 

= Y.'^T ^Y.'^xGT , (1) 

1=2 1=2 

with n = 7V-|- 1, vj the interaction between the projectile and the I target subsystem, and G is the propagator 

G = (e+ - K -Y,Vxjy\ (2) 

Here, E is the total energy and is related to the total incident kinetic energy Ei = 2^in\ by = -Ei + eo where 
eo is the target ground state energy. At this stage we use non-relativistic kinematics. The inclusion of relativistic 
kinematics will be discussed later. The operator K corresponds to the total kinetic energy of the projectile and N 
target sub-systems in the projectile-target center of mass frame. In Eq. Q Vjj is the interaction between subsystems 
X and J. Equivalently, we may write the propagator in terms of the kinetic energy operator for the projectile Ki in 
the center of mass of the interacting projectile-target (P-T) system, and the target nucleus Hamiltonian, ffo 

G = {E+ - Ki - Hq)'\ (3) 

The total transition amplitude, Eq.Q, can be rewritten as 

T^Y.^x + Y.^tGY,tj + --- . (4) 

where the projectile-X subsystem transition amplitude tj is given by 

Ti = vi + vxGtj . (5) 

We note that the propagator in tj contains the target Hamiltonian and thus it is still at this stage a many-body 
operator. In the limit when the target nucleus subsystems are weakly bound to each other, the multiple scattering 
expansion to the P-T transition amplitude is expected to converge rapidly 0, 0| and the MST expansion Eq. Q can 
be used. 

We apply next this formalism to proton scattering from a target of two and three bound structureless subsystems. 

III. THE SINGLE SCATTERING 3-BODY PROBLEM 

We first consider the scattering of a projectile of mass mi from a target (such as ^^Be) assumed to be well described 
by a two body model with two subsystems (of valence particle and a core) labeled here as 2 and 3 of masses m2 
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and ma respectively. Let ki,k2,k^, (k'-^,k2,k'^) be the initial (final) momenta of the projectile and the two cluster 
subsystems. 

Neglecting core excitation, the target wave function is 

*(^,6) = ['/'23(r)®<,93(6)] , (6) 

where (/33(f3) is the core internal wave function and 023 is the wave function describing the relative motion of the 
(2,3) pair. 

The elastic transition amplitude to first order in the projectile-subsystem transition amplitudes is 

T = (fcl$|r2|fci$) + (fci$|r3|fci$) + --- (7) 

There are two approaches in handhng the dynamics of the few-body system: The impulse (I) and the Fixed Scatterer 
(or frozen halo, or adiabatic) (FS) approach. They are related in high energy regime for special cases of the scattering 
amplitudes. 



A. The factorized impulse approximation [FIA] 



Within the impulse approximation, the interaction between the clusters Vxj is assumed to have a negligible dy- 
namical effect on the scattering of the projectile from the individual target subsystems and therefore can be neglected. 
The operator projectile-I target subsystem transition amplitude rj is then replaced by 



ti = vj + vxGotj 
where Go contains only the kinetic energy operator K 



(8) 



(9) 



The transition amplitude is still a many body operator, because the kinetic energy operator has contributions from 
the projectile and all Af target subsystems. The interaction between the target subsystems leads to terms in 3rd 
order of the tx transition amplitude. As we shall see, this projectile-X subsystem amplitude can be reduced, after 
suitable approximations, to a free two-body amplitude evaluated at the appropriate energy. Accepting the validity of 
replacing n by tj in Eq. I^J we obtain the multiple scattering expansion 



T 



lA 



(10) 



In the single scattering approximation (SA) only the first term is taken into account. 

Let us consider the scattering from subsystem 1 = 2. Because from the dynamical point of view we want to reduce 
the problem to the projectile scattering from each subsystem, we take as relevant coordinates the relative momentum 
between projectile and subsystem X=2, (ji^2, and the relative momentum between the target subsystems, g2,3, defined 
here as 



91,2 = 



7712 fci — mifc2 



M- 



12 



92,3 



7713^2 - rn2kz 



M- 



23 



(11) 



where M12 = nii + 1712, M23 = m2 + m^. In the three-particle cm. frame, Vt = X]i=i ki = Q {V't = X]i=i ^'i — 0); 
and the propagator Go is 

-1 



Go — 



E+ 



E^ 



2fll2 
2^12 



112 



912 



2Atl2,3 
2^*12,3 



923 + 71— fci 

JW23 



H -1 



(12) 



with ^12 = mim2/Mi2 and /ii2,3 = m3Mi2/A/i23, Afi23 = mi+m,2 + m^. The total energy, neglecting binding effects, 
is 



E = 



2fi 



1,23 



-kl 



2/Z: 



-Q 



23 



2 

23 5 



(13) 
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with /ii,23 = W1M23/M123, and (523 an average relative momentum between the target subsystems, to be specified 
below. The single scattering matrix elements are 

(fci023|i2|fcl023> = J dq23(f>23 (^'?23 " 'J^^^ 

X (91,21*2(0)12)191,2)023 (923) • (14) 
Here, A is the projectile momentum transfer 

A = k[-ki . (15) 
The initial and final relative momenta between the projectile and the subsystem 1 = 2, qi,2 and q[ 2 are respectively 

qi,2 = ki 92,3 = P12K1 - ^1292,3 , 

Ml,23 "^2 

9l,2 = «1 92,3 = Pl2fc 1 - "1292,3 

Ml, 23 "^2 

/^i2 Aii2 ^ , M12 "l3 X 
-«i 92,3 + 71—^ 



Ml, 23 W2 ' m2 M23 
= /3i2fci - Q!i292,3 + ai27i2A , (16) 

with /3i2 = ^12/^1,23, ai2 = M12/TO2, and 712 = /i23/w2. We shall use these whenever a simplified notation is required. 
These new parameters satisfy P12 + ai27i2 = 1, from which it follows, together with Eq. II16|I . that the condition 
q'l 2 — 9*1,2 = A necessarily holds. The energy parameter i2!i2 is 

^2 / ^ N 2 

(1,12=^; 9l3 + ^fci ■ (17) 

Ml2,3 \ -'^23 / 

The single-scattering matrix elements of Eq. Ijl4|l involve a full folding integral of a product of a transition amplitude 
and a target form factor. This integral may be quite involved. As we shall see below, in the high energy regime, 
the relative momentum of the interacting pair 9*23 can be approximated by a suitable value Q23 • Once replaced in 
Eqs. II16I17|I one obtains a factorized impulse approximation (FIA) expression. 



T''^{E) = (Q'i2|i2(c^l2)|Ql2)p23 (^^A 

+ (Q'l3</'3|i3(^13)|'/'3Ql3)p23 (^^A^ , (18) 

where Qu is approximate relative momentum projectile-subsystem I, ujn the energy parameter and p23 the target 
form factor 

P23 (Ai) = J ^923023 (923 " 023 (9*23) (19) 

In the Hmit of a heavy subsystem ma ^ 7712, we obtain the expected Hmit of P23 (i^^) — 1- Equivalently, we 

can write the scattering amplitudes in terms of the transferred momentum K12 = Q 12 ^ Q12 and total momentum 
^12 = [Q'i2 + Qi2]/2: 

T™(£;) = i2(tc'12,Kl2,-^12)p23 

+ t3(wi3,Kl3,-€l3)p23 f^A^) . (20) 



The scattering amplitudes are on-shell if \Qii\ = IQ'ul = y^^^ijoju/h. In this case the scattering amplitudes are 

local. As we shall show below, although 9i^2~9i,2 = A, the approximate FIA may not satisfy K12 — A. The scattering 
amplitudes are related to the transition amplitude according to 

\FiE)\' ^ ^^ki^msr (21) 
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FIG. 1: (Color online) Kinematics for the 3-body scattering in the impulse approximation of Kujawski and Lambert (KL). 
The dashed lines represent the momentum transfer for the case of the 3-body system (dashed dark line) and the 2-body 
system (dashed light line). The vectors fc(fc') represent the initial (final) projectile momenta, and Q{Q') the initial (final) 
projectile-subsystem T relative momenta. 



for projectile target scattering, and 



hvi 



IX 



|ix(wii)|' 



(22) 



for projectile subsystem I scattering. The elastic differential angular distribution is then evaluated from the total 
scattering amplitude on-shell 



A 



23 



'23 



( ^K] , 
VM23 ' 



(23) 



where the normalization factors are 



IT 



1 



diJii /dE 



(24) 



We now discuss several factorized impulse approximations that can be found in the literature. All these approximations 
use Eq. ifnji as starting point. The differences between the models arise from the approximations performed to obtain 
a factorized expression of the form II18|I . 



1. The 3b on-shell FIA [Kujawski and Lambert] 

In the factorized on-shell approximation, discussed in the work of Kujawski and Lambert (KL) jlo|, the incident 
and outgoing relative momenta of the subsystems 923(923) ^re neglected by comparison with the projectile momentum. 
That is, inserting 

IQ23I = (25) 
into Eq. Hlfill7|l . one obtains for the relative momenta between the projectile and subsystem 1 = 2 

Q,^ = J^k, , Q'^^ = J^k'^ , (26) 
Ml, 23 Ml, 23 



and for the energy parameter 

^'^^ ~ I I 

2a<12 M-1,23 



C.12 = ^\Qi2\' = ^E . (27) 
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Similar expressions are obtained for the scattering from subsystem X = 3 



and for the energy parameter 



Qi3 = i^fci , Q[^ = J^k[ , (28) 

Ml, 23 Ml, 23 



ui3 = -^E . (29) 

Ml, 23 



The kinematics for the KL approximation is represented schematically in Fig. Clearly, in this approach the relative 
projectile-subsystem J momentum is a fraction of the total transfered momentum A 



Kii = -^A . (30) 

Ml, 23 



The single scattering terms are then 



Ml, 23 Ml, 23 \-lW23 

+ iJ^k[\Uco,s)\^ki)P23(^A) . (31) 

Ml, 23 Ml, 23 \M23 J 

By construction from Eqs. Ij2(il27ll . the matrix elements of the transition amplitudes t2 and ^3 are on-shell. Equivalently 
one may write the transition amplitudes as a function of the projectile momentum transfer 

„ r / M12 x\ ( r 
T = ^2(^^12, A)p23 Ti— A 

Ml, 23 \A'-/23 



+ ^3(^^13, A)p23 77—^ 

Ml, 23 V-/W23 



(32) 



Using Eqs. l|3Tll and lf23jl . the total scattering amplitude F can be written in terms of the scattering amplitudes for 
each subsystem /j, with normaHzation coefficients 

Mi= , 1=1,2 . (33) 

V Mil / 



2. The 3b on-shell FIA [RihanJ 

The optimal factorized approximation discussed in the work of Rihan 0, ^| was formulated for a three body 
problem in the context of the multiple scattering expansion of the optical potential. It was to our knowledge never 
applied to a specific scattering problem. In this approximation, the relative momentum between the two subsystems 
5*23 is taken as the mid-point where the product of the wave functions peaks in Eq. Ijl4|l . that is. 



1 rns 



123 



2M- 



23 



(34) 



Substituting this into Eq. ifTfiji . the relative momenta are 



Q12 = /3i2fci 



"12712 



O' /Q U' < "12712 r 

Q12 = P12K1 H 7, — A 



(35) 



and the energy parameter is 



t^l2 



2mi 



IS 



12 



Ml, 23 



(36) 
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2-body 




FIG. 2: (Color online) Kinematics for the 3-body scattering in optimal factorization impulse approximation of Rihan. The 
lines and vectors have the same meaning as in Fig. Q 



with 



Cl2 



A*l,23 
. Ml2 . 



[I + C12+C12 COS 6*] 



(37) 



Here, 9 is the projectile-target scattering angle and 

C12 = -ai27i2 + ^"12712 • (38) 

For small scattering angles C12 ~ 1. Similar expressions can be obtained for the scattering from subsystem X = 3. 

The kinematics for the Rihan approximation is represented schematically in Fig. [2l Within this approximation, 
using /3i2 + q;i27i2 = 1, it follows that 

Kii = A (39) 

The single-scattering matrix elements are 



T = (/3i2fci H 7, — A|t2(t^i2)|/3i2fci — A)p23 71— A 



2 " \M: 



i23 



+ (/3l3fcl + ^^A|£3(-13)|/?13fcl-^A),23(^Aj . (40) 

Equivalently one may write the transition amplitude in terms of the projectile momentum transfer 



T = h (-12, a) P23 (^a) h (c.13, a) P23 (^a) . (41) 

By construction from Eq. Ij35|l and Eq. Ij36|l . the matrix elements are on-shell. 

Making use of Eqs. ll3Tll and ll23|l . the total scattering amplitude F can be written in terms of the scattering 
amplitudes for each subsystem /x, where the normaHzation coefficients are 

Ml, 23 \ / 1 



We note that for small scattering angles the energy parameters reduce to Eqs. ll29ll and ll27|l and the normalization 
coefficients to Eq. l(33jl . 
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FIG. 3: (Color online) Kinematics for the 3-body scattering in the on-shell impulse approximation of Chew. The lines and 
vectors have the same meaning as in Fig. Q 



3. The 3b on-shell FIA [Chew] 

For completeness we discuss now the approximation discussed in [T^ , although no practical application will be 
made here. Within this approximation, the relative incident momenta between the subsystems 923 is again neglected 
when compared with the incident projectile momentum. That is, setting (523 = one gets 

Qi2 = -^fci . (43) 

A'1,23 

In addition, it is also assumed that 

|Q'i2l = |Qi2| , l«l = |A|. (44) 

The kinematics in the Chew approximation is represented schematically in Fig. El 

The energy parameters, 1012 and cjia are given by Eq. Ij27|l and Eq. H29II respectively. By construction from these 
equations the matrix elements are on the energy shell. The normalization coefficients for the angular scattering 
amplitudes are identical to those derived in the KL approximation Eq. l(33jl . 



4- The 3b on-shell FIA [Crespo and Johnson] 

In this approach, followed in the most current applications of MST the initial internal relative momentum 

g23 is neglected in the transition matrix elements, and thus Eq. is satisfied. In addition, 

Q'i2 = Qi2 + A (45) 

The energy parameters, W12 and wia are given by Eq. lf27jl and Eq. lf29jl respectively. 

The single-scattering terms can then be written as a function of the transition amplitude for proton-subsystem 
scattering and the density for the subsystem P23, 

r=(^fcl+A|t2(^12)|-^^l)p23 (^^] 
Ml, 23 Ml, 23 V^W23 / 

Ml, 23 Ml, 23 V-'^^23 / 

(46) 

or equivalently 

r==i2(wi2,A,/Ci2,(/))M23 {j^^^ 
+ 4(W13, A,/Ci3,0)M23 (^-^A^ , 

(47) 
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2-body 




FIG. 4: (Color online) Kinematics for the 3-body scattering in the off-energy shell, OFFES, impulse approximation. The lines 
and vectors have the same meaning as in Fig. Q 



with 4) the angle between A and /C. On the energy shell, (j) — 7r/2. The matrix elements of the transition amplitude 
are half on the energy shell since Q12 = \/2/ii2Wi2/?i, but Q'12 7^ Q12, as represented schematicaly in Fig. 31 

This approach is however impractical if one has only on-shell scattering amplitudes for the scattering from the 
subsystems. This is the case, for example, when these amplitudes are obtained by fitting scattering data. We note 
that for small scattering angles, the scattering angle for projectile subsystem I ~2 scattering, 6*12, satisfies 

Ail(23) „ 



n,2 



(48) 



This means that the range of physical momentum transfers over which the on-shell ti2 transition amplitude is defined 
is smaller than the range of momentum transfers accessible in for the scattering from the 3-body system. 

If the dependence of the transition amplitude tx{ujix, A,/Cii, 0) in Eq. (j47ll on the total momentum is assumed to 
be small then we can replace the matrix elements by their on-shell values at the appropriate energy. On the energy 
shell, the total momentum JC is related to the momentum transfer A, and (j) = 7r/2. In most current multiple scattering 
applications 0,13,0 the on-shell approximation is made within a partial wave expansion of the transition amplitude, 
but the angular direction is kept constant. We refer to this factorized impulse approximation approach as projected 
on-energy shell [POES]. 

Making use of Eq. (|17jl and Eq. lf23jl the total scattering amplitude F can be written in terms of the scattering 
amplitudes for each subsystem fx, where the normaHzation coefficients are given as in Eq. 



B. The Fixed scatterer or adiabatic approximation [FSA] 

Within the fixed scatterer or adiabatic approximation, the internal Hamiltonian between the clusters is taken to a 
constant H, that is, the operator projectile-I target subsystem transition amplitude tx is replaced by 

ix = vi + vxGoix (49) 

where Go 

Go = {E+ - Ki - . (50) 

Even after this simplification, the solution of this problem involves the solution of a system of coupled equations, 
which in practice requires truncation of the angular momenta in order to make the problem solvable fl^l • A notable 
simplification of the problem is achieved in general when one neglects the interaction of the projectile with all the 
fragments except the core 0,^^, the so called core- recoil model. In this particular case, the few-body problem can 
be exactly solved, leading to a simple factorized expression, the product of the transition amplitude for the scattering 
from the core evaluated at the appropriate energy and a structure form factor. 

Within the FSA framework for the case where the interaction between the projectile and fragment 3 is neglected 
the total transition amplitude takes the form, 

T^^^ = {k[\i,{E)\k,}p,,(^^^A^ . (51) 
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We note that in the FSA approach, the energy parameter is E = ?iA;^/2/ii 23 and thus, a distinctive feature of the FSA 
is that the two-body ampUtude is calculated with the reduced mass /Ui_23, instead of the reduced mass appropriate 
for the two-body scattering which appears in the FIA scattering amplitude. In the general case, i.e., when the 
interaction with the two target subsystems is considered, one can perform a multiple scattering expansion of the 
few-body amplitude in terms of the individual T-matrices tx for the fragments which, in leading order, yields the 
factorized expression, 

+ {k[\i,iE)\h)p,,(^^^A^ , (52) 

In principle, the fixed scatterer approximation is conceptually different from the factorized impulse approximation, in 
any of the versions discussed above. However, we shall show further down that in the high energy limit, the calculated 
elastic scattering observables using the FSA and the FIA (FOES or Rihan) are essentially identical. The formal 
relation between the FSA and FIA and their relation with the Glauber theory will be explored elsewhere. 



IV. THE SINGLE SCATTERING 4-BODY PROBLEM 



We consider now the case of proton scattering from a nucleus assumed to be well described by a three body model 
(for example ^^Li), as a core (here labeled as subsystem 4) and two valence weakly bound systems (subsystems 2 and 
3). Neglecting core excitation, the target wave function is 

R, U) = ['/'23,4(r, R) ^4(^4)] , (53) 

where (p4{£,4) is the core internal wave function and 4>23,4{r,R) the three body valence wave function relative to the 
core. The internal degrees of freedom of the core are denoted globally as ^4. The total transition amplitude is then 

T'^ = (fc;$|£i2|fcl$) + (fci$|tl3|fcl$) + (fcl$|tl4|fcl$) 

+ ••• (54) 

where within the single scattering approximation only the first three terms arc taken into account. We first consider 
the situation where particle 1 scatters from one of the valence systems, named subsystem 2. The relevant Jacobian 
coordinates are 

m2ki-mik2 

^^■^ = W2 ^''^ 

7723^2-^2^3 

'''' = M23 ^''^ 

.-23,4 = , (57) 

JW234 

where P12 = ki + ^2, P34 = ^3 + ^4, etc. In the center of mass of the total four-body system, Pt = Xli=i — Tfie 
intermediate states propagator is given by 

/ fi2 t2 *2 \-i 

Go = (e^ - - - ^9?2,34 ) (58) 

\ 2/(Xl2 2/i34 ^Ml2,34 / 



with 

M34P12 - M12P34 

912,34 = (59) 

J"1234 

7774^3 - 7773 ^4 

'^^■^ = W4 

and /7i2, /734, /Ui2,34 = M12M34/M1234 the appropriate reduced masses. Using the fact that the operator ti2 is 
independent of the spatial variables of the core, the single scattering matrix elements are given by 
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(^'l 023,41^2 |fcl</>23,4) 



^9*23 ^923,4023,4 9*23 



7713 



A, 923,' 



m4 



^ I (9 1 2l*2(ti'12)|gi, 2)023,4 ((f23, 923,4) (61) 



M23 ' M234 

The initial (and final) relative momentum between the projectile and the subsystem 1 = 2, (fi 2 (9 '1.2) a-re respectively 

TO2 



91,2 



9 '1,2 



M12 




M12 


923 + - 

1 


Ml, 234 


7712 


M12 




M12 


9 23 + 


Ml, 234 


?7l2 


M12 


h' 


M12 


9*23 + 


Ml, 234 ^ 


m2 


M12 


TO3 


m2 TO4 


m2 


.M23 


M23 M234 



-923,4 



23 
TO2 



M23 

A , 



923,4 



(62) 



with ai2 = M12/W2 defined as in the 3-body case. The energy parameter (Di2 is 



0^12 — E - 



2m 



34 



TO3M234 
M23M; 



34 



-923,4 



7714 



1 2 



34 



-923 



2/i 



12,34 



M34 r , ™2 ^ , ^ ' 

Kl + T:;— 923,4 + (723 



M- 



234 



M2: 



(63) 



We note that expressions Eq. Ij62|l and Eq. Ij6^^|l reduce to the 3-body problem in the limit where we take nii = 0. 

We now consider the situation where particle 1 scatters from the core named subsystem 4. The relevant Jacobian 
coordinates are 92,3,9*23,4 and 



91,4 



Al- 



ii 



The intermediate states propagator is given by 



Go 



2fi 



-<ll4 



with 



14 



914,23 



2m: 



23 



-9I3 



2^ 



14,23 



9 14.23 



M23P14 - Mi4Pi 



23 



M- 



1234 



(64) 



(65) 



(66) 



and Mi4, M23; Mi4,23 — M14M23/M1234 the appropriate reduced masses. The single scattering matrix elements for the 
scattering from the core 



(fcl02304|i4|<i54023,4fcl) = / / ^9*23 <^923, 4023,4 ( 9*23,9*23,4 + JJ^^ ) (91,404 1*4(^14) 19549*1,4)023,4 (9*23 , 9*23,4)(67) 



^23, 923,4 -1- , , 
V ^#234 

The initial (and final) relative momentum between the projectile and the subsystem 1 = 4, 9i_4 (9 4) are respectively 



91,4 



9*1.4 



Ml4 ^ , Ml4 ^ 

ki H 923.4 

Ml. 234 "^4 



Ml4 



'^11 9 23.4 • 

Ml, 234 TO4 



The energy parameter 0)14 is 



LUl4 = E — 



2pL 



23 



-92^3 



2^ 



14,23 



23 



kl - 9*23,4 



234 



(68) 



(69) 
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In the limit where we take = 0, these equations reduce to the 3-body case. As in the 3-body case, it is desirable 
to obtain after suitable approximations a factorized expression of the transition amplitude matrix elements and the 
target form factor 

X 023,4(Ql,Q2) . (70) 

In here 023.4 (Qi, Q2) is the Fourier transform of the wave function of the two body valence system relative to the core 
023^ (r, R). The factorized impulse approximation expression for the 4-body case is 

T^^^ - (Q'l2|i2K2)|Ql2>P23.4f^A,-^A 



VM23 ' M: 



234 

(Q'l3l4K3)|Ql3)p23.4 f^A,-^A 
\i\i23 -'W234 

(Q'l404|4(wi4)|</'4Q14)P23,4 fo, ^A 

\ JW234 



(71) 



where Q and uj are the appropriate relative momenta and energy parameter respectively to be discussed bellow, which 
are a generalization of the 3 body case. 



1. The 4b on-shell FIA [Kujawski and Lambert] 

Within this approximation, the relative momenta between the subsystems 923 and 923,4 are neglected when compared 
with the incident projectile momenta in the matrix elements of the projectile subsystem transition matrix elements 
whenever they appear, that is, 

IQ23| = |Q23l=0 

IQ23,4| = |Qi3,4l=0 . (72) 

From Eq. l(fi2|l one obtains for the relative momenta between the projectile and subsystem 1 = 2 

A M12 ^1 M12 f, 

0.12 = Kl , Ql2 = fci (73) 

Ml, 234 Ml, 234 

and for the energy parameter 



wi2 = (74) 

Ml, 234 



and similarly for the scattering from subsystems X = 3,4. The single scattering terms can then be written as 

T = (^fc;|£4(c.u)|^fcl)p23,4 [ 0, ^a) 
Ml, 234 Ml, 234 \ jV-'234 / 

+ (^A:i|t2(c.,2)|^fci)M23,4 f^A,-^A) 

Ml, 234 Ml, 234 \JW23 -("234 / 

+ (^A:;|t3(.,3)l^fcl>M23,4f^A,-^A) . 
Ml, 234 Ml, 234 \-'W23 ^"234 / 



(75) 



By construction, from Eqs. Ij 731 7411 the matrix elements of the transition amplitudes are on-shell. 

As in the three-body case the total scattering amplitude, -F, can be written in terms of the scattering amplitudes 
for each subsystem /x, where the normalization coefficients are 

Mx=f^V ,1 = 2,3,4 (76) 
V Mil / 
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2. The 4b on-shell FIA [RihanJ 

The extension of the optimal factorized approximation discussed in the work of Rihan (Til IT^ to the four-body 
problem is straightforward. The relative momentum between the subsystems is taken to the mid-point value where 
the product of the wave function peaks. For the scattering from subsystem 1 = 2 this yields 

0.3 = 



Q 



2M23 
1 7714 



23,4 



2M, 



A 



234 



(77) 



which leads to 



^ A 7^ "12712 

Ql,2 = Puki 7^ — A , 



(78) 



with ai2 defined as in the three-body case and (3i2 = ^^12/^1, 234, 712 — A/34/M234. Note that these expressions reduce 
to the 3-body case in the limit where we take TO4 = 0. Since this approximation involves on-shell matrix elements the 
energy parameter can be evaluated as 



(^12 



2^1 



12 



Ml, 234 



(79) 



with 



C12 



Ml, 234 
fJ.12 



[1 +C12 +C12 C0S( 



(80) 



where 6 is the scattering angle and 



C12 = -ai27i2 + 2"i27i2 



(81) 



The normalization factor is formally identical to the three-body case, Eq. II42II . with C12 replaced by Eq. Ij8()|l . Similar 
expressions can be derived for the case of the scattering from subsystem 1 = 3. In the case of the scattering from 
subsystem 1 = 4, the optimal approximation prescribes 



Q23 = 

V23,4 — ^ 



2M' 



which leads to 



234 



^ Ml4 ^ Ml4 1 A/23 X 
2l.4 = fcl XTT A 

Ml, 234 m4 2 M234 
ai47i4 



(82) 



-A 



(83) 



with /3i4 



714 = Af23/A^234 and for the energy parameter. 



Wl4 



2/i 



14 



eJ^C44 
Ml, 234 



(84) 



with 



Ml, 234 
Ml4 



Ci4 cos ( 



(85) 



where 



Cl4 



-ai47i4 -I- -a?47i4 



(86) 
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3. The 4b on-shell FIA [Crespo and Johnson] 



In this approach the initial relative momenta (723,923,4 are neglected in the transition matrix elements momentum 
transfer. For the scattering of subsystem T — 2 

Ql2^-^k , Q'l2 = Sl2 + A (87) 
Ml, 234 

and similarly for the scattering of the other subsystems. 



V. THE "He STRUCTURE MODEL 



The ^He is here described as a three-body system n+n+'^B.e. The bound wave functions are obtained by solving the 
Schrodinger equation in hyperspherical coordinates. We consider two structure models defined in terms of different 
effective 3-body (3B) potentials, which are introduced to overcome the underbinding caused by the other closed 
channels, most important of which the t+t breakup. In both models the n-^He potential is taken from Ref. [77,1^, 
and use the GPT NN potential 0| with spin-orbit and tensor components. In the first model (R5) of 0^ the 3B 
effective potential in the hyperspherical coordinates is given as a function of the hyperradius p 

In the second model (R2) of jlJl the potential is defined as 

V^'^'Xp) = -U3cM-p') ■ (89) 

In these equations p — p/po, where po= 5 fm and 1 fm for R5 and R2 respectively. The strength of the 3B effective 
potential is tuned to reproduce the experimental three-body separation energy, with V3=-1.60 MeV and U3=-293.5 
MeV. The models R5 and R2 predict, with an a particle rms matter radius of 1.49 fm, ^He rms matter radii of 2.50 
fm and 2.35 fm respectively. 



VI. RESULTS 



In this section we evaluate the elastic scattering differential cross section for the scattering of protons on ^He at 
Eiah = 717 MeV within the four-body FIA single scattering approximation and Fixed scatterer approximation FSA. 
In the case under study, the single scattering terms involves contributions from the valence nucleons and from the 
'^He core. 



A. p-|-cluster 



In the factorized impulse approximations discussed in the present work, the total T-matrix is expressed in terms of 
the on-shell matrix elements of the two-body amplitudes evaluated at an appropriate momentum transfer and energy. 
The NN on-shell scattering amplitudes were obtained from a realistic NN Paris interaction, as in 0,0, 01 • 

The two-body scattering amplitude for the scattering p+'^He was calculated with a phenomenological optical po- 
tential, with parameters obtained by fitting existing data for the elastic scattering p+'^He at Ep = 700 MeV 2^ 
and 800 MeV 2 1|. We consider two different parametrizations for the p + a optical potential. In the first one, based 
on the work of Baldini-Neto et al. we take the optical potential as a sum of two terms, real and imaginary, of 
Gaussian shape 

U{r) = V^Tir, ro) + iW^T{r, n) (90) 

with T{r, Tx) = exp(— r^/r^). The values of the radii and the real and imaginary depths were adjusted simultaneously 
in order to minimize the with the experimental data. These fits were performed with the computer code FRESCO 
[2^, version frxy. The final values of these parameters where Vb=+73 MeV, ro=1.37 fm, Wo=-162 MeV and ri=1.32 
fm. 



15 




FIG. 5: (Color online) Calculated elastic scattering for p+'^He at £"^=700 MeV (solid lines) and 800 MeV (dashed lines), using 
a Gaussian and a Woods-Saxon parametrization as described in the text. The data are taken from Refs. [2I. l2^. l24l| 



We have also taken a phenomenological potential of Woods-Saxon (WS) shape, i.e.: 

U{r) = Vofir,ro,ao) + tWJ{r,n,a,) (91) 

with f{r,r^,ax) = [1 + exp[(r - i?^)/a^]~\ where Rx ^r^A^I'^. 

To reduce the number of degrees of freedom in the best fit search we constrained the parameters to ro=ri and ao=ai. 
With this condition, a best-fit analysis of the data yields the values Vb=+72 MeV, Wi=-115 MeV, ro=ri=1.14 fm, and 
ao=ai=0.31 fm. Note that with both geometries the real potential was positive in our fits, indicating a dominance of 
the repulsive part in the p+'^He interaction at these energies. The data and calculated differential elastic cross sections 
obtained with these parametrized potentials for p+^He at Ep=700 and 800 MeV are shown in Fig. [3 The solid and 
dashed lines, corresponding to i?p=700 MeV and 800 MeV respectively, reproduce very well the forward scattering 
data for both the Gaussian and Woods Saxon parametrization. By contrast, the calculated differential cross section 
tends to deviate from the data of [l^l at larger angles in the case of the Gaussian parametrization. The WS optical 
model improves the fit in this angular region. 

B. p+^He 

The two-body NN and p + a scattering amplitudes obtained by fitting elastic data are now used to evaluate the 
total transition amplitude for p+^He scattering. 
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-t (Gey/cf 

FIG. 6: (Color online) Calculated p+^He elastic scattering at 717 MeV for two different optical potentials for p+'^He. The data 
is taken from [27ll2SH2flll. 



The Coulomb interaction was included in an approximate way as summarized in appendix A. The elastic scattering 
observables shown in here were evaluated using relativistic kinematics as discussed in detail in appendix B. The 
relativistic kinematic effects nevertheless were found to be small. In Fig. El we study the dependence of the calculated 
differential cross section for p+^He elastic scattering with respect to the underlying p+'^He optical potential as a 
function of the squared four momentum transfer —t = A^. The solid and dashed lines represent the FIA-POES 
calculation using the Gaussian and WS Optical model parametrizations for core scattering. At higher momentum 
transfers the WS calculation gives a slightly bigger reduction of the cross section. The difference with respect to the 
Gaussian parameterization is small, and both parametrizations predict essentially the same differential cross section. 
The conclusions of the present work are therefore essentially independent of the underlying OM potential for the 
scattering from the core. We shall be using for definiteness the WS potential in all subsequent calculations. 

In Fig. 13 the thick solid line represents the calculated differential elastic cross section for p+^He using the FIA- 
POES approximation. The dashed Hue was obtained neglecting the single scattering contribution from the valence 
neutrons. By comparing these two calculations, one finds that the core contribution dominates the large angle region 
but underpredicts the forward differential cross section. Also shown in the figure is the calculated elastic cross section 
for p+^He (thin solid line). For values of -t < 0.03 (GeV/c)^ {9^.^. < 11°) the calculated p+^Re cross section is 
significantly bigger than the p+^He cross section, providing a good agreement with the experimental data. Comparison 
of the two FIA calculations indicate that this enhancement on the cross section is mostly due to the presence of the 
valence neutrons, since the FIA calculation with core contribution alone is very close to the p+'*He distribution at 
these angles. 

By contrast, for larger values of momentum transfer, the few-body calculation is smaller than the p-f^He curve. One 
notes from the figure that the two FIA calculations (core and core+valence single scattering) are very close to each 
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FIG. 7: (Color online) Experimental and calculated p+ He elastic scattering at 717 MeV within the FIA-POES approximation. 
For comparison purposes, the OM calculation for p+''He elastic scattering, at the same energy per nucleon, is also included. 



other at large momentum transfer, indicating that, at least within the single scattering approximation, the depletion 
of the cross section at larger angles is mostly a consequence of core recoil effects. Although the experimental p+^He 
data exhibit a depletion of the cross section with respect to the p+^He distribution, the reduction predicted by the 
FIA calculation is too small to explain the data. 

In Fig.lHlwe compare the calculated p+^He cross section using the FIA-POES (thick soHd Hne), the FIA-KL (dashed- 
dotted line) and the FSA (adiabatic) approximation (long dashed). The cross section calculated with the FIA-Rihan 
approximation is very similar with the FIA-POES and therefore it is omitted from the figure. The difference between 
the predicted differential cross sections largely refiects the different rms matter radii. 

The calculated differential cross section using FIA (POES and Rihan) and FSA are very similar. The elastic 
scattering observable calculated using the KL impulse approximation decays more slowly and the agreement with the 
data is comparatively even poor at the large angle region. 

We now study the sensitivity of the results with respect to the ^He wavefunction. In particular, we compare in 
Fig- El the FIA-POES calculations obtained with the R5 (thick soHd fine) and R2 (dashed Hne) ^He models discussed 
in sectionEl Both models give nearly identical results at low momentum transfer, but the R5, having a larger radius 
than R2, exhibits a larger reduction of the cross section at large scattering angles. However, these differences do not 
appear to explain the disagreement between the calculations and the data, and are much smaller than the differences 
obtained when using different scattering approximations. 

From the analysis presented in this section, we can conclude that the factorized impulse approximations, in the 
different versions here presented, describe fairly well the elastic scattering data for p+^He at small momentum transfer. 
In particular, they all show an increase of the cross section with respect to the p+'^He scattering, showing the effect 
of the valence neutrons on the elastic scattering. At momentum transfers above 0.05 (GeV/c)^ all the calculations 
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FIG. 8: (Color online) Experimental and theoretical p+^He elastic scattering at 717 MeV per nucleon, for several scattering 
models discussed in the text. The calculated p+*He elastic distribution is included for comparison. 

tend to overestimate the data, irrespective of the structure model used. 

VII. CONCLUSIONS 

In this work we have reviewed and compared several approaches for the scattering of a nucleus by a weakly bound 
composite, based on the multiple scattering expansion of the total scattering amplitude, namely, Factorized Impulse 
and Fixed Scatterer/Adiabatic Approximations. 

The Factorized Impulse Approximation (FIA) neglects the inter-cluster interaction and approximates the relative 
momentum between the loosely bound clusters. As for the Fixed Scatterer approximation (FSA) it replaces the 
internal Hamiltonian by a constant. 

As a common feature, all these approaches express the single-scattering term as a product of a two-body scattering 
amplitude for the scattering of the fragment times a structure formfactor, which depends on the internal wavefunction 
of the composite system. This factorized form provides a simple interpretation of the scattering observable, by 
separating the role of the structure from the reaction dynamics. 

These approaches have been then appHed to the scattering of ^He on protons at E=717 MeV per nucleon, using 
p-n amplitudes derived from a realistic NN potential, and p+'^He amplitudes obtained from an optical model fit of 
existing elastic data. All the approximations here succeed in reproducing the p+^He elastic data at small scattering 
angles. We found it crucial to include the effect of the valence neutrons to the single scattering contributions in order 
to explain the increase of the cross sections with respect to the p+'^Ue elastic data at the same energy per nucleon. 
At larger angles, both FSA and FIA calculations tend to overestimate the data. 
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FIG. 9: (Color online) Experimental and calculated p+^He elastic scattering at 717 MeV for two different models for the ®He 
ground state. 

This disagreement at larger angles could be traced to the neglect of higher order terms in the T-matrix expansion 
and will be discussed elsewhere. 
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APPENDIX A: TREATMENT OF THE COULOMB POTENTIAL 

The elastic differential angular distribution is evaluated from the total scattering amplitude 

\M23 -'W234 



+ <73K...,).,„(£A,-g-A) 

1/2 ^ / -^^23 
+ /4(t^l4,Ki4)p23,4 0, — A 

V JW234 



(Al) 

In this equation, the scattering amplitude from the core includes the Coulomb interaction. The central scattering 
amplitude is taken as 

pi I 



fm = fcio) 



-^^exp(2za,.)[(i + l)Tf+(Ar) 

W14 — 



L 

+ LTI-{N)]Pl {cos 9) . (A2) 

In Eq. ljA2p . f^iO) is the Coulomb scattering amplitude due to a point charge Z^e and Qu is the asymptotic 
projectile-core wave number. The T^^{N) are defined according to 

Gxp [2i(5f=^(7V)] - 1 
T^^{N) ^ , (A3) 

where L± denotes the orbital and total angular momenta, J = i ± | and 5^^{N) are the Coulomb modified phase 
shifts It follows from Eqs. (|AllA2ll that the point Coulomb scattering amplitudes appears multiplied by the 

structure form factor P23,4 'M^^ • Although this is an approximate treatment of the Coulomb interaction, it will 
only have effects at very small projectile momentum transfer and therefore does not modify the conclusions of the 
present work. 

APPENDIX B: RELATIVISTIC KINEMATIC EFFECTS 

In this section we describe how the expressions of sections III and IV should be modified in order to take into account 
relativistic kinematics. For simplification we consider the 3-body case Eq. (II8II . Expressions can be straightforward 
generalized to the 4-body case. For simplicity, we shall take in this section c = 1. 

Within MST, the projectile-target scattering scattering amplitude is constructed from the projectile-subsystem 
scattering amplitude as described in the text. We shall discuss the relativistic kinematics for the scattering of both 
the composite and each subsystem. 

Let us consider the elastic scattering of projectile labelled 1 with a target A with rest mass mi and ttia respectively, 

1 + A^l + A (Bl) 

In our case A represents the composite (2+3) two-body system. In relativistic kinematics one introduces the Lorentz 
invariant Mandelstam variable sia = [Ei + EaY ill- The differential cross section for projectile-target with respect 
to the four momentum transfer tiA is related to the cm. differential cross section as: 

da t: da 

di^^^dn ' ^ ' 

where tiA — — |Ap — 2fcJ(cos0c.m. — 1) and fci the projectile momentum in the projectile-target CM frame, that can 
be obtained from the Mandelstam variable sia- The cross section is evaluated from the scattering amplitude, related 
to the matrix elements of the total transition amplitude according to P| 



21 



with w = y/s\A — TOi — 'n^A- 

Let us now consider now the projectile scattering from subsystem 1=2. As before, one introduces the Mandelstam 
variable for the interacting pair ,si2. For simplification, we shall take the relative momentum of the interacting pair 
q'23 nonrelativistically. We shall consider in here the two situations where 523 = (as in KL and POES impulse 
approximation) and 523 7^ (as in optimal approximation discussed by Rihan) 

• In the former case, in the laboratory frame A;2 = ^3 = and the Mandelstam variable S12 can be readily 
evaluated 

si2=ml + ml + 2m2 (rm + T^^^^) . (B4) 

• In the later case, where 923 = it is more convenient to evaluate S12 in the projectile-subsystem 1 = 2 
CM frame. In this case, 



S12 



(^^Jml+kl + ^m^2 + kl^ -{k + kf (B5) 



To evaluate S12, one takes k\ from the Mandelstam invariant s\a- The last term gives 

(fcl + k^f = \ll2kl [1 + COS0e.m.] (B6) 



and 



k'2 — k\ 



X12 + ^7i2 - 712X12 cos6ic.r 



(B7) 



where X12 = m^jMii + 5712 



The differential cross section for the scattering from subsystem X with respect to the four momentum transfer tyi , is 
related to the cm. differential cross section as: 

dh2 Qf2dfli2 ' ^ ' 

where ti2 = — |Ki2p = 2Qf2(cos0i,2 — 1) and Q12 the projectile momentum in the projectile-subsystem CM frame 
that can be obtained from S12. The elastic scattering amplitude is related to the T-matrix elements through: 

/2 = ^fc?^(Q;2|i2(u;i2)|Ql2) , (B9) 

with a;i2 = \/si2 — mi — m2 ■ 



[1] C. J. Joachain. Quantum collision theory. North-Holland, 1987. 

[2] E. O. Alt, P. Grassberger, and W. Sandhas. Phys. Rev., C 1:85, 1970. 

[3] A. C. Fonseca. Phys. Rev., C 30:35, 1984. 

[4] T. Matsumoto, E. Hiyama, K. Ogata, Y. Iseri, M. Kamimura, S. Chiba, and M. Yahiro. Phys. Rev., C 70:061601, 2004. 

[5] M. L. Goldberger and K. M. Watson. Collision Theory. John Wiley and Sons, New York, 1964. 

[6] K. M. Watson. Phys. Rev., 105:1338, 1957. 

[7] R. Crespo and R. C. Johnson. Phys. Rev. 60, 034007, 1999. 

[8] R. Crespo, I. J. Thompson, and A. A. Korshcninnikov. Phys. Rev., C 66, 2002. 

[9] R. Crespo. Proceedings of rnb5, 3-8 april 2000, divonne, france. Nucl. Phys., A701:429c, 2002. 
[10] E. Kujawski and E. Lambert. Ann. Phys., 81:591, 1973. 
[11] T. H. Rihan. Phys. Rev., D 1235:1235, 1977. 
[12] T. H. Rihan. Phys. Rev., C 53:2328, 1996. 
[13] C. F. Chew. Phys. Rev., 84:1057-1058, 1951. 
[14] G. F. Chew. Phys. Rev., 80:196, 1950. 

[15] J. A. Christley, J.S. Al-Khalili, J. A. Tostcvin, and R. C. Johnson. Nucl. Phys., A 624:275, 1997. 
[16] R. C. Johnson, J. S. Al-Khalili, and J. A. Tostevin. Phys. Rev. Lett, 79:2771, 1997. 
[17] J. Bang and C. Gignoux. Nucl. Phys., A 313:119, 1979. 



22 



[18] I. J. Thompson, B. V. Danilin, V. D. Efros, J. S. Vaagen, J. M. Bang, and M. V. Zhukov. Phys. Rev., C 61:24318, 2000. 
[19] P. Pires D. Gogny and R. de Tourreil. Phys. Lett, 32B:591, 1970. 

[20] B. V. Danilin, I. J. Thompson, M. V. Zhukov, and J. S. Vaagen. Nud. Phys., A632:383, 1998. 
[21] N. K. Timofeyuk and I. ,J. Thompson. Phys. Rev., C 61:044608, 2000. 
[22] G. D. Alkhazov et al. Phys. Rev. Lett, 78:2313, 1997. 

[23] O. G. Grebenjuk, A. V. Khanadeev, G. A. Korolev, S. I. Manayenkov, J. Saudinos, G. N. Velichko, and A. A. Vorobyov. 

Nud. Phys., A 500:637, 1989. 
[24] L. G. Arnold, B. C. Clark, and R. L. Mercer. Phys. Rev., C 19:917, 1979. 

[25] E. Baldini-Neto, B. V. Carlson, R. A. Rego, and M. S. Hussein. Nud. Phys., A 724:345, 2003. 
[26] I. J. Thompson. Comp. Phys. Rep. 7, 167, 1988. 

[27] F. Aksouh. Ph.D. Thesis. Investigation of the core-halo structure of the neutron rich nuclei ^'^ He by intermediate- energy 

elastic proton scattering at high momentum transfer. PhD thesis, 2002. 
[28] P. Egelhof. Invited talk presented at the International Symposium on Physics of Unstable Nuclei. Nud. Phys., A 722:C254- 

C260, 2002. 

[29] F. Aksough et al. Proceedings of the 10th International Conference on nuclear reaction mechanisms, Varenna, 2003. Review 

of the University of Milano, Ricerca Scientifica ed educazione permanente, Suppl. 122, 2003. 
[30] R. Crespo and J. A. Tostevin. Phys. Rev. C 41, 2615, 1990. 



